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We elaborate a modification of the deformable two-sublattice Mitsui model of [Levitskii R.R. et al., Phys. Rev. B. 
2003, 67, 174112] and [Levitskii R.R. et al., Condons. Matter Phys., 2005, 8, 881] that consistently takes into 
account diagonal components of the strain tensor, arising either due to external pressures or due to thermal 
expansion. We calculate the related to those strains thermal, piezoelectric, and elastic characteristics of the 
system. Using the developed fitting procedure, a set of the model parameters is found for the case of Rochelle 
salt crystals, providing a satisfactory agreement with the available experimental data for the hydrostatic and 
uniaxial pressure dependences of the Curie temperatures, temperature dependences of spontaneous diago- 
nal strains, linear thermal expansion coefficients, elastic constants eg and c^, piezoelectric coefficients du 
and gii {i = 1, 2, 3). The hydrostatic pressure variation of dielectric permittivity is described using a derived 
expression for the permittivity of a partially clamped crystal. The dipole moments and the asymmetry param- 
eter of Rochelle salt are found to increase with hydrostatic pressure. 
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1. Introduction 

The Mitsui model [1] (two-sublattice Ising model with asymmetric potentials) has been origi- 
nally proposed for description of the reentrant phase behavior in Rochelle salt crystals. It considers 
the motion of certain ordering units in two interpenetrating sublattices of asymmetric double-well 
■otentials. The model with certain modifications is also applicable to several other ferroelectrics 
like those of the Rochelle salt type (deuterated and ammonium-doped Q Rochelle salt at least 
at low doping ), RbHS04 type |3,i|, AgNa(N02)2 (SSN [i,|3|), SASD type (NaNH4S04-2H20 
and NaNH4Se04-2H20), etc. 

Very often, inclusion of deformational effects into the Mitsui model is indispensable for a proper 
description of the system behavior even at ambient pressure. Thus, a strong piezoelectricity asso- 
ciated with polarization Pi and shear strain £4 essentially affects the dynamic dielectric response 
of Rochelle salt due to the effect of crystal clamping by the high-frequency measuring field. The 
conventional Mitsui model yields a qualitatively incorrect behavior of the relaxation time and 
dynamic dielectric permittivity near the Curie temperatures. This problem is resolved [9|] by tak- 
ing into account the piezoelectric coupling with £4, also permitting to describe the phenomena of 
piezoelectric resonance and sound attenuation pO*!. 

The sequence of phase transitions observed in NH4HSO4 crystals can be described within 
the mean field approximation for the Mitsui model only with temperature dependent interaction 
constants [5|. It means that one must take into account the effect of thermal expansion, which is, 
ultimately, a deformational effect. 

Since high pressure studies are the only means to continuously vary the system geometrical 
parameters, such as the orientation angles of atomic groups that form the dipole moments, inter- 
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atomic distances, hydrogen bonds parameters, etc, as well as the interparticle interactions, and 
other parameters of the system, they can provide a valuable information on the mechanism of the 
phase transitions in ferroelectric crystals. A better insight is obtained if the effects of hydrostatic 
and various uniaxial and biaxial pressures are explored. Irrespectively of the crystal symmetry, 
these pressures produce diagonal components of the lattice strain tensor Si {i = 1,2,3). In low- 
symmetry systems, shear strains Sj {j — 4, 5, 6) can be induced as well. The diagonal strains also 
arise due to the thermal expansion of the crystals. 

Nowadays. Rochelle salt ^11,-, 14] and other crystals [31 described by the Mitsui model often serve 
as test materials for experiments with nanosize phenomena. Properties of the ferroelectric nano- 
inclusions in a solid matrix are strongly affected by surface tension and thermal mismatch stresses. 
Thus, radial stresses in the plane perpendicular to spontaneous polarization arise in Rochelle salt 
nanorods grown in pores of alumina films [ill [l^ , inducing diagonal strains only. A phenomeno- 
logical theory of ferroelectric properties of such nanorods was presented in [l^ . 

The simplest way to incorporate an external pressure into a spin model is to consider its param- 
eters (interaction constants) phenomenologically as linear functions of pressure. This approach has 
been used for different versions of the Mitsui model to describe the hydrostatic pressure variation 
of the transition temperatures in Rochelle salt [17] and SASD [18]. However, to describe a uniaxial 
stress dependence of Tq in the same way, one would have to find new values of the fitting param- 
eters for each stress direction. A unified model description of hydrostatic and uniaxial pressure 
effects, and, at the same time, the crystal thermal expansion at ambient pressure, with a single set 
of the fitting parameters is possible if we include into the model the lattice diagonal strains, instead 
of the pressures. Such a microscopic-like model of bulk crystals that includes the diagonal strains 
will be very helpful when one develops a model description of the above mentioned nanocrystal 
behavior. 

The goal of the present paper is to develop a modification of the deformable Mitsui model with 
the shear strain £4 [9|, which would also take into account the diagonal components of the lattice 
strain tensor. The first attempt to create such a modification was made in (3]. The interaction 
constants and the asymmetry parameter were taken to be linear functions of the diagonal strains. 
Expressions for the piezoelectric and elastic characteristics of Rochelle salt, associated with these 
strains, have been obtained. However, all actual calculations were performed in the approximation 
of zero thermal strains; the external hydrostatic or uniaxial pressure effect was not considered, and 
the fitting procedure was inappropriate. Here we shall use the model of [l9| and take into account 
the thermal expansion strains properly. We shall develop a consistent fitting procedure, free from 
the drawbacks of the previous [19j calculations, allowing us to obtain the above mentioned unified 
description of high pressure effects, thermal expansion, as well as physical properties of Rochelle 
salt associated with the diagonal strains. 

Bulk Rochelle salt undergoes two second-order phase transitions at Tci — 255 K and Tc2 = 
297 K, with the intermediate ferroelectric phase. Spontaneous polarization Pi is directed along the 
a axis, accompanied by spontaneous shear strain £4 in the he plane. The crystal is orthorhombic 
(space group P2i2i2i) in the paraelectric phases and monoclinic (P2ill) in the ferroelectric phase. 
As it follows from the analysis of symmetry elements of its point group 222 and of the uniaxial 
pressure point group oo/mmrn, no uniaxial or biaxial pressure applied along the orthorhombic 
crystallographic axes of a Rochelle salt crystal changes its symmetry. Neither do the hydrostatic 
pressure or thermal expansion of the crystal. 

The mechanism of phase transitions in Rochelle salt remains rather obscure. According to the 



most recent measurements |20|, |2l| , the largest displacements at a ferroelectric phase transition are 



undergone by the Os, O9, Oio oxygens. It appears that the dipoles of the Mitsui model, moving in 
the double-well potentials, should plausibly be associated with the OH9 and OHio groups. Their 
motion, coupled with displacive vibrations of OHg groups seems to be responsible for the phase 
transitions, as well as for the spontaneous polarization. It would be interesting to elucidate the 
pressure variation of the dipole moments, as well as of the asymmetry parameters of the double- well 
potentials. This can shed some light on the details of the transition mechanism. 

The set of experimental data for Rochelle salt, used for verification of the present modification 
of the Mitsui model, includes the hydrostatic pressure dependences of transition temperatures 
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(2^ . [23I and dielectric permittivity f23| , the uniaxial 2^, 2^ and biaxial j2d\ pressure dependences 
of the Curie temperatures. Related to the diagonal strains, components of piezoelectric (e.g. du, 
i — 1, 2, 3) and elastic (0^4) tensors ji^-H^i appearing in the ferroelectric phase only, should be also 
described by the model. The other characteristics, included into the fitting, are thermal expansion 
coefHcients and dilatations [.30,,^], as well as the diagonal-strain-related elastic constantsjlH Cij 
or compliances Sij — 1,2,3). 

We introduce the diagonal strains ei, £2, £3 into the two-sublattice Mitsui model with the shear 
strain £4 P| in the manner it was done in [l^. In section 2 we make a further modification of the 
model by taking into account the host lattice contributions into the thermal expansion. Section 3 
contains the obtained expressions for thermal, elastic, dielectric, and piezoelectric characteristics. 
In section 4 we propose a new fitting procedure. Using the found set of the fitting parameters for 
the Rochelle salt crystals, we show that the developed theory is capable of describing the entire 
complex of the phenomena, related to the diagonal strains: thermal expansion, temperature behav- 
ior of monoclinic piezoelectric and elastic characteristics, hydrostatic, uniaxial, and biaxial pressure 
dependences of the Curie temperatures in Rochelle salt. The behavior of the dielectric permittiv- 
ity under hydrostatic pressure is described using the derived expression for the permittivity of a 
partially clamped crystal. 

2. System thermodynamics in the presence of diagonal strains 

We consider an orthorhombic piezoelectric crystal in the paraelectric phase, to which an external 
hydrostatic or uniaxial or biaxial (along the crystallographic axes) pressure can be applied. All these 
pressures produce diagonal strains £i (i = 1, 2, 3); these strains are also present due to the thermal 
expansion. Paraelectric piezoelectricity is associated with the shear strain £4. The Rochelle salt 
symmetry is presumed having the spontaneous polarization directed along the a axis and coupled 
to the strain £4. 

We start with the modified two-sublattice Mitsui model with a piezoelectric coupling to the 
shear strain £4 and with the diagonal strains [l^ . A unit cell of the model consists of two dipoles 
(two sublattices), oppositely oriented along the a-axis (the axis of spontaneous polarization); they 
are compensated in the paraelectric phases and get non-compensated in the ferroelectric phase. 
The actual unit cell of a real Rochelle salt crystal is twice as large and contains four dipoles. At 
the transition to the monoclinic phase the angle between the b and c axes changes from 7r/2 to 
7r/2 — £4. The diagonal strains £1, £2, £3 describe the relative changes in the lattice constants a, b, 
and c, respectively, due to thermal expansion or under pressure. 

In the mean field approximation, the model Hamiltonian reads 



9 



where 77/ = {cqf); N is the number of the unit cells; J, K are the Fourier-transforms (at k = 0) of 
the constants of interaction between pseudospins belonging to the same and to different sublattices, 
respectively. 

The phenomenological part of the Hamiltonian [/scod is a "seed" energy of the host lattice of 
heavy ions which forms the asymmetric double- well potentials for the pseudospins. For the case of 
Rochelle salt symmetry in presence of diagonal strains, shear strain £4, and field Ei, it reads 

3 3 

(^scod - yC44 £4 - Voe^^Eihi ^Xll^l + y % ^i^3 ~ ^'0 2^ % Qfj (J ~ 1, )£-, . (2) 

ij"=l ij=l 

Here £o is the vacuum permittivity; vq is the unit cell volume of the model. The three first terms 
in Useed are the elastic, piezoelectric, and electric contributions due to the shear strain £4 and 
longitudinal electric field Ei. Two last terms are related to the diagonal strains. Here cf^, c^°, ef^^ 
are the "seed" constants describing the phenomenological contributions of the crystal lattice into 
the corresponding observed quantities c^, cfj, and 614. In fact, the index E in cf^*^ is redundant, as 
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the difference between the observed cfj and cfj at i,j = 1, 2, 3 is neghgible. The "seed" quantities are 
zeros if the corresponding observed quantities are zeros in the most symmetric phase (orthorhombic 
in the case of Rochelle salt). 

The last term, absent in the earlier model [l^, is the contribution of the host lattice into 
the energy of thermal expansion, are the "seed" thermal expansion coefficients; T° are the 
temperatures at which the components of this contribution vanish. It is known that the thermal 
strains can be set to be equal to zero at any arbitrary temperature Tq (the reference point for 
thermal expansion). This can be achieved by choosing the values of accordingly; they will 
differ from Tq due to the pseudospin system contributions to the thermal expansion. To take into 
account the "seed" contribution of the host lattice is indispensable for a proper description of 
thermal expansion. 

The coefficients 

^^(1) ^ \jm + \Km + A - 2i/;4£4 + E{2) = i Jry2 + ^i^r/i - A - 2^4e4 + Mi^^i (3) 

in (HI) are the local mean fields acting on pseudospins of the first and second sublattices in the q^^ 
unit cell. The parameter A describes the asymmetry of the double well potential; is the effective 
dipole moment. The model parameter ^/;4 describes the internal field created by the piezoelectric 
coupling with £4 and essentially determines the piezoelectric and elastic characteristics associated 
with the shear strain £4 [1, [l^. It is also assumed that a longitudinal electric field Ei is applied. 

3 

J±K = Jo±Ko + 2Y,4^fe^. (4) 



as well as the asymmetry parameter 



3 



Ao + X! ^3z£» (5) 

i=l 



are taken to be linear functions of the diagonal strains [l9'|. Here -0^ are introduced simply as 
the expansion coefficients. However, for J and K such an expansion is equivalent to taking into 
account the electrostrictive coupling with the diagonal strains. For A it implicitly describes the 
changes in the asymmetry parameter due to the changes produced by external pressure or thermal 
expansion in the interatomic distances and in the geometric parameters of the potential, like the 
distance between the potential wells, etc. The parameters 7/13^, are analogous to the deformational 
potentials 27 — dA/de introduced in the Anderson- Halperin-Varma-Phillips [SS, 34| model of two- 



level systems with asymmetric double-well potentials in order to describe ultrasound attenuation 
and thermal conductivity in amorphous solids at very low temperatures. 

The thermodynamic potential of the considered model is obtained within the mean field ap- 
proximation in the following form 



r rr.. rr J + ^ ^2 J^K^ 2 lu 2 1, , 7 + ,1^5 A 

92E{j)i,T) ^Usccd^ — £, H — cr In cosh cosh hwo^^Piei, (6) 

where /3 — l/k^T, is the Boltzmann constant, and 

7-/3(^^-|^e-204£4 + mi?i), S = 13 (^^^-^a + A 

Here pi — P2 — Ps — Ph for hydrostatic pressure; pi ^ and P2 ~ Pa = for the uniaxial pressure 
applied along the axis a, etc. The shear pressure p4 is introduced formally, in order to find the 
elastic and piezoelectric characteristics associated with it; after that it is put equal to zero. 
We introduced the following linear combinations of the mean pseudospin values 
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^ is the parameter of ferroelectric ordering in the system. The parameters ^ and a are determined 
from the saddle point of the thermodynamic potential ([B]): a minimum of g2E with respect to ^ 
and a maximum with respect to cr are realized at equilibrium. The corresponding equations are 

sinh 7 sinh 5 



^ cosh 7 + cosh 5 ' cosh 7 + cosh 5 

3. Physical characteristics of Rochelle salt related to diagonal strains 



(7) 



Using the following thermodynamic relations 



1 



'2E 



= 0, 



El 



2E 



dEi 



where v = vo{l + X]i=i ^0 pressure and temperature dependent unit cell volume, and 

retaining only linear in Si terms in the "seed" contributions, we obtain expressions for strains and 
polarization 



£4 = 



^0.+a^(T-Tf) + ^^.f (V+C2 + V^7a2 + 2V.3,^), (* = 1^3), (8) 
2^-4 



P4_ 

,E0 
L-44 



^1 



Pi = e°4£4 + eoXn^^i + — ^ ■ 

V 



(9) 
(10) 



Here s,^° are the elements of the matrix inverse to the matrix of "seed" elastic constants cf"P. 

Two first terms in equation ([5]) represent the host system contributions into the Hooke's law 
and thermal expansion with regular pressure and temperature behavior. The sum in equation ([5]) 
gives the pseudospin subsystem contributions into the strains, having anomalous behavior in the 
ferroelectric phase. The second term in equation ([SJ was absent in the previous model [19] . 

From equations ([5|)- (fTU)) we can derive expressions for other characteristics related to the diago- 
nal strains. Thus, the coefficients of linear thermal expansion are obtained in a rather cumbersome 
form 



1 



fc=i 



>^ E 4° [^^^(^2^ - ^37) + i^'-T + ^3,)(A27 - ^5<5)] , 



(11) 



where 



B = 

Cjk ~ ~ 
<^5 = Ai 



/3 



2vq[lp2 ~ A(^3) 



A 



J is a unit matrix. The other notations are 



(^2 = 1- 

Ai = 1 - 



^Ai-/3 ^^(Ai 



A2 - 2^a, 



A-- 



(P3^ Xi + 



-A - ^2j: 



ml 
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Alternatively, the coefRcients of thermal expansion can be found by numerical differentiation of 
equation ([5]) for the strains with respect to temperature; the results, of course, coincide. The 
coefficients ai are expected to have small anomalies at the Curie temperatures (sol . 3f |. 

The molar specific heat at constant pressure of the model is obtained from the molar entropy 



NA/ ^glE{e^,T) 



2 V 



dT 



VqNa 



ij=l 



R 
2" 



2 In 2 + In cosh cosh 7^ — (Jcr 



where giE{si, T) = g2E{'^i, T) —vq '^^Pief. R is the universal gas constant, and A'a is the Avogadro 
constant. Thus, the molar specific heat of the model is 



vqNaT 



dT 
(A27 - (fsS) 



ij=l 
3 



R 



4((/32 - A(^3) 



i=l 



(12) 



As we shall see, it has small anomalies at the transition points. To obtain the total specific heat 
of a crystal that can be compared to experimental data, we have to add to equation p2l) a regular 
term linear in temperature (within the considered temperature range) that would correspond to a 
contribution of lattice vibrations not taken into account within our model. Thus, 



Cp Cvibr 1 



Cvibr = A + BT; 



(13) 



the coefficients A and B will be specified by fitting equation (IT51) to experimental data. Often, an 
inverse procedure is performed, when a regular linear contribution is subtracted from the experi- 
mental data; the obtained result is then compared with the theoretical specific heat of the ordering 
subsystem. 

The other found characteristics are, in particular, the elastic constants at constant electric field 
= 1,2,3) 



S4 



\dej / E,T 

EO _ P 

2«o^2 

^^4 iPii 



as well as the monoclinic piezoelectric coefficients 



(14) 
(15) 



dPi 
( d£i 



Ei,T 



\dEi 



V \ 2ip2 



J2 



i + Eti£. 



(16) 



where sfl is the matrix of elastic compliances, inverse to the matrix of elastic constants c^, and 

I3{K + J) 



= Ai - 

The other piezoelectric and elastic characteristics are 



fdEi 



Pi 



£0X1 



gii 



fdEi 



(A? - A^) 



du 



£0X1 



dpi 



+ eiihij 
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Here, 



is the dielectric susceptibility of a clamped crystal, and 

3 



eo\dEi) 2veo ip2 - A(p3 eq 



-Vei.dii (18) 



i=l 



is the static dielectric susceptibility of a mechanically free crystal [19|. Here we introduce the 
following notations 

fj-i - Hi - ^?/'4ai4 , «i4 - -go , Xii - Xii H ei4C(;^4 . 

C44 £0 

In paraelectric phases, this expression for the free su scep tibility coincides with that obtained within 
the modified Mitsui model without thermal strains [10|. The sum e^^ ^ Gudu , different from zero 
in the ferroelectric phase and not exceeding 5% of the total susceptibility, was absent in the earlier 
model. 

As one can easily verify, monoclinic quantities en, du, hu, gu, cf^, cf^ [i — 1,2,3) differ from 
zero only at non-zero polarization, in agreement with the symmetry considerations. 

The temperature of the second order phase transition is determined from the condition that 
dielectric susceptibility of a free crystal Xii diverges at T — Tc. From equation (llSp and using 
equation ([7]) we obtain 

where the model parameters J. K, A are taken at Tq, being renormalized by diagonal strains 
according to equation Q. 

Equation (fT9|) is valid both for the ambient pressure case and for the stressed crystal. It can be 
rewritten in the two following convenient forms: 



\ 



1 (20) 

4 vncf° 



which gives an explicit expression for a at the transition points, and 



} ^ Os.Sci + CTc V'i £c^ = - Aq cTc + 2fcBTcArccosh4 



4fcBrc vocf^kBTc 
useful in the fitting procedure. Here Ed are the strains at the Curie temperature. 



4. Numerical calculations 

4.1. Fitting procedure 

The model parameters must provide a fit of the theory to the experimental data for the following 
characteristics: the Curie temperatures at ambient pressure Tck (fc = 1,2 in Rochelle salt) and 
their hydrostatic and uniaxial pressure slopes dTck/dph and dTck/dpj, the temperature curves 
of thermal expansion strains Ei, linear thermal expansion coefficients, monoclinic piezomodules 
gii, and elastic constants Cij and ca = 1,3). Simultaneously we check for an agreement with 
experiment for the quantities previously described [l3 by the modified Mitsui model without 
thermal strains, such as spontaneous polarization, static free and clamped dielectric susceptibilities 
Xii : piezomodule (ii4, specific heat, elastic constant at constant field c|^, as well as microwave 
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dielectric permittivity £ii{v, T). A detailed analysis of the effect of diagonal strains on the physical 
characteristics of Rochelle salt associated with the shear strain £4 will be given elsewhere. 

The adopted values of the model parameters are given in table [1] details of the fitting procedure 
are described below. 



Table 1. The model parameters used for description of Rochelle salt. 



do 


/in 




(K) 


A,, Jh-rr, 




(10-^" C-m) 




Ci4 

C/m^ 


All 


0.3162 


0.662 


764.63 


1476.46 


745.14 


-750 


8.7 


-0.0008 


0.033 


10.1 








7/)j~/fcB 


(K) 




1p3l/kB 


V'32/fcB 


i'ss/kB 




-13000 


-12000 


-9850 


10614 


13537 


1080 


-12125 


-15993 


-6043 




a? 


(10-^ K-^) 




Tl 


(K) 


rpO 










5.800 


3.353 


4.333 


238.35 


185.32 


311.71 










Cll 


C12 


Cl3 


(10'° N/m 






C44 








2.842 


1.794 


1.541 


4.219 


2.031 


3.987 


1.18 









The temperature variation of the order parameter ^ is determined by numerical minimization 
of the thermodynamic potential ([5]); cr is found from equation ([7]); the strains are determined from 
equations ([S]) and As the reference point for thermal expansion (where Si — at pi — 0) we 
chose the upper transition temperature Tc2 = 297 K. This condition allows us to express from 
equation ([5]) via a^, ip~ , and ip^i. Strictly speaking, are not the fitting parameters of the model, 
since the reference point can be chosen arbitrarily. At 308 K and at ambient pressure, the lattice 
constants are ^ a = 11.927 A, 6 = 14.292 A, c = 6.225 A. 

The "seed" linear thermal expansion coefficients were chosen to yield the thermal strains 
Ei at Tci equal to af'^^{Tci — 7c2), where af^^ are the experimental [31| values of the expansion 
coefficients in the middle of the ferroelectric phase (at 275 K). 

Special care has been taken that below 20 kbar for the hydrostatic pressure and below 200 bar 
for uniaxial or biaxial pressures and between and 350 K, no additional phase transition takes 
place in the system, apart from those taking place at ambient pressure. 

The number and (if any) temperature and order of the phase transitions for the Mitsui model 
without thermal strains are usually analyzed in terms of the dimensionless variables a and b 



K - J ^ 8A 

"''4^44 



K + J+^^sfr ^ K + J+^^^l^^^ ^^^^ 



VqT I'i VQ 

and the dimensionless transition temperature 



t-o^ 4fcBrc (22) 



The phase diagram of the conventional (undeformable) Mitsui model in the (a, b) plane |36l - l38| 
shows the regions with different numbers and types of the phase transitions; its topology is not 
changed by inclusion of the shear strain £4. It has been found that only in a very narrow region 
of the (a, b) plane, the system undergoes two second order phase transitions with the intermediate 
ferroelectric phase. 

In the presence of diagonal strains, a and b become functions of temperature and pressure. In 
the fitting procedure we shall deal with the values of a and b at the upper Curie temperature and 
at ambient pressure do and bo- Absence of additional phase transitions at the chosen values of 
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the model parameters is verified directly, by calculating the order parameter at all temperatures 
between and 350 K and at pressures below 25 kbar (hydrostatic) or 200 bar (uniaxial). The values 
of ao and 69 should be from the same region of the (a, b) phase diagram of the undeformable Mitsui 
model that yield two second order phase transitions. With decreasing do, the maximal values of 
spontaneous polarization, spontaneous strain £4, and anomalous parts of diagonal strains increase. 
We choose the value of ao that gives the best agreement with experiment for these characteristics. 
Once the values of ao, bo, tp4, and cf^ are chosen, we are in a position to find Jo, Kq, and Ao, 
using equations (gHl), (ED), and (P^ . 

The parameters cf^, ip4, and /ii are varied around their values obtained in the previous study 
[o'l in order to get the best fit for spontaneous polarization Pi, piezoelectric coefficient c?i4, static 
free and clamped Xi'i dielectric susceptibilities and dynamic dielectric permittivity eii{v,T). The 
dipole moment fii is assumed to decrease linearly with an increasing temperature as 

[i + kTiT-Tc2)]. 

The values of ^5 and fcx are given in table [1] 

We require that the theoretical values of the elastic constants (i,j = 1,3) at Tc2 should 
coincide with their experimental values |3S] , available for 307 K (this is a reasonable approximation 
due to a very weak temperature dependence of c^). Thus, we can easily determine cf^^ using 
equation (fH)) . 

It is required that the best possible description of experimental [2^ . [25| uniaxial pressures 
dependence of the two Curie temperatures should be obtained. To this end, at the chosen ao, 
bo, tp4, C44 , , the six parameters and ip^i are determined from six linear equations (j2ip 
written at T^j, = T^j, + {dTck/dpijPi {k = 1,2 and i = 1,2, 3) at pi = 100 bar and combined with 
equation l|20[). The slopes dTck/dpi were varied around their experimental values [23|. The strains 
at the Curie temperatures were approximated as 

during the fitting. The obtained values of 6^ and were found to be independent of the used 
values of pi . 

One of parameters, say, ip'^ , can be determined from the condition that the lower Curie 
temperature should be Tci = 255 K. For the two remaining parameters tp^ and ipt there is the 
condition that the two calculated transition temperatures at hydrostatic pressure of 1 kbar and 
the lower transition temperature at 20 kbar would be in agreement with the experimental data 



22l. I23j. However, the dependences of the transition temperatures on uniaxial pressures and on 
the hydrostatic pressure are not completely independent, (this will be discussed later). Therefore, 
normally ipi and can be varied continuously in certain ranges and provide correct theoretical 
dependences Tck{Ph) at given ao, 60 and V'4- From these ranges we should select the values which 
yield the best fit for the temperature curves of the piezoelectric coefficients gu, anomalous parts 
of thermal strains Ssi in the ferroelectric phase, and thermal expansion coefficients a^. It should 
be mentioned, however, that a perfect fit both for 171^ ([l^l) and Ssi ([Hj) cannot be obtained 
simultaneously, so a certain compromise has to be made. 

A criterion to make an unambiguous choice of the theory parameters hardly exists. Since the 
chosen set of ipf^ is not unique, it is not possible to precisely establish the temperature and pressure 
variation of the interaction constants. However, the overall tendency is such that the hydrostatic 
compression enhances the asymmetry parameter A, as well as the constants of interactions between 
the pseudospins within the same and in different sublattices. Pressure slopes of A. J, and K are 
very sensitive to the choice of bo at given do, while the observable quantities are not that much 
sensitive. The average slopes are about 1-3%/kbar for J and 4.5-7%/kbar for A and K. This issue 
will be explored in more detail elsewhere. 

4.2. Thermal expansion and specific heat 

The temperature dependence of the diagonal strains caused by thermal expansion of a crystal 
in the absence of external pressures is plotted in the inset to figure [TJ 
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The experimental points, obtained from the data for thermal dilatations [3l|, are well described 
by the proposed theory. 




— I ' 1 1 1 1 1 ' 1 ^iu -I 1 ' 1 ' 1 1 1 1 1 

260 280 300 320 340 240 260 280 300 320 



Figure 1. Anomalous parts of thermal diago- 
nal strains of Rochelle salt as functions of tem- 
perature: 1, □: £si; 2, O- £s2; 3, A: £53. Inset: 
total strains as functions of temperature: 1,0: 
£1; 2, O- £2; 3, A: £3 Lines: a theory; symbols: 
experimental points taken from [Slj. 



Figure 2. Coefficients of linear thermal expan- 
sion of Rochelle salt as functions of tempera- 
ture. 1, □, ■: ai; 2, Q, •: 0:2; 3, A, A: 0:3. 
Lines: a theory. Open and closed symbols are 
experimental points taken from [3l| and from 
|40| . respectively. 



In the ferroelectric phase, the ei{T) curves have small bucklings (anomalous parts) caused by 
electrostrictive coupling to spontaneous polarization. In order to extract these anomalous parts of 
the diagonal strains, Imai [31| extrapolated the measured temperature curves of the strains from 
the paraelectric phases onto the ferroelectric phase and subtracted them from the total measured 
strains. With the same purpose, we calculate some hypothetical paraelectric strains (coinciding 
in the paraelectric phases with the actual strains Si) from equation ([5]) by putting ^ = at all 
temperatures and determining cr from equation ([7]), and find the spontaneous strains as Esi = £i—£^. 
The obtained results are shown in the major part of figure [TJ As one can see, at the adopted values 
of the model parameters, the theory well reproduces the asymmetric shape of the £s3(T) curve, 
but underestimates the magnitude of Sgi and £52- The agreement can be improved by choosing 
different values of the model parameters, but the agreement with gii{T) will be spoiled. 

The corresponding linear thermal expansion coefficients are shown in figure [5] The theoretical 
and experimental values of their jumps at the Curie temperatures are summarized in table [2j 



Table 2. The calculated jumps of the thermal expansion coefficients and specific heat at Curie 
temperatures. The values in parentheses are experimental data of [sij. 





Aa2 
(lO-^K-i) 




Acp 
(J/mol K) 


Tci -6.0 (-7.0) 
Tc2 4.7 (7.1) 


3.0 (3.5) 
-2.1 (-3.7) 


3.7 (4.5) 
-1.1 (-2.0) 


0.74 
1.09 



We have a fairly good agreement with experiment [3l| for ai and 03, both for the signs and for 
the values of the coefficient anomalies at the transition points, as well as for the temperature slope 
of the curves in the ferroelectric phase. The agreement is worse for a2- It should be noted that 
the experimental behavior of a2(T) is somewhat different from that of ai and 013, with a marked 
decrease below the lower Curie temperature and a large difference between the coefficient values 
in the two paraelectric phases. No perceptible temperature variation in the ferroelectric phase was 
experimentally detected, also in contrast with the ai and behavior. The theoretical temperature 
curve of a2(T) is very much like that of and qualitatively similar to that of ai. The theoretical 
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a2 in the upper paraelectric phase weU accords with the single experimental value of [40'| obtained 
from the synchrotron radiation Renninger scan, although the reported error of these measurements 
is so large that the values of [sj] also fall in this error range (see figure [2]). 



Figure 3. Specific heat of Rochelle salt as a function of temperature. Lines: a theory. Solid line: 
this work; dashed line: the modified Mitsui model without the thermal strains [9|]. Symbols: 
experimental points taken from 43]. 

Figure [3] shows that the present model yields a better agreement with experimental data for 
the small anomalies of specific heat of Rochelle salt at the Curie points than it was obtained with 
the earlier model [qI , especially for the magnitude of the upper anomaly. The regular contribution 
of lattice vibrations was taken to be Cvibr = 105.845 + 0.855T (J/mol K) for the present model and 
103.456 + 0.944r (J/mol K) for the previous model. 

The calculated values of the specific heat jumps are given in tablej^l We obtain positive anoma- 
lies at both transition points, in accordance with the most recent measurements ^43j . The jumps 
were defined as the differences between the ferroelectric and paraelectric values of specific heat or 
expansion coefficients at the Curie points. 

4.3. Diagonal-strain-related piezoelectric and elastic constants 

Figure m shows temperature dependences of piezoelectric constants gu. As one can see, a fairly 
good agreement with experiment is obtained, including the opposite signs of gu and of gi2 and 
^13, asymmetric shape of (713 (T) dependence, as well as the interception of the 1712 and 1713 curves 
near the lower Curie point. The overall behavior of gii{T) is similar to that of esi{T). 
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Figure 4. Piezoelectric constants of Ro- 
chelle salt as functions of temperature at 
Si = 5 kV/cm: 1, □: <?ii; 2, Q: 312; 3, A: gig. 
Lines: a theory; symbols: experimental points 
taken from [2711. 



Figure 5. Elastic constants C24, + and -f 
of Rochelle salt as functions of temperature 
at E\ = 500 V/cm. Lines: a thcOTv; symbols: 
experimental points taken from [29| . 
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We do not depict the calculated elastic constants cfj and cfj {i, j = 1, 2, 3), as they are practically 
temperature independent between 230 and 330 K. A very small variation can be detected in the 
ferroelectric phase for cfl, with the difference between cfl and constant cf, being less than 1% of 
at most. The maximal values of monoclinic constants cf^ are more than by one order of magnitude 
smaller than cfj . The shape of the experimental + vs T curve is qualitatively reproduced by 
the proposed model, as seen in figure [Sj A quantitative agreement with experiment is reasonable. 

The coefficients of piezoelectric strain dn are shown in figure IH] at zero and high bias fields 
El. The data of [27] were extracted from the given therein values of the dugn product and of 
gii. In absence of external field du actually diverges at the Curie temperatures, due to the term 
proportional to sf^ in equation (jl6p ; the bias field smears out these anomalies and lowers the peaks 
of the coefficients. A good agreement with the experimental points is obtained. 

£/j.(10"''C/N) 
200 

100 





-100 



-200 



250 275 300 325 350 375 

Figure 6. Piezoelectric coefRcients du of Rochelle salt as functions of temperature. The inset: 
dii as a function of temperature at different fields Ei (kV/cm): 1, ; 2, A: 0.304; 3, ■: 5. Lines: 
a theory; symbols: experimental points taken from (27| (■) and [j^ l (A). 



d„(10'-C/N) 




4.4. Hydrostatic pressure effects 

Below we shall discuss how the high-pressure effects are described by the proposed modification 
of the Mitsui model. Figure [7] shows the calculated hydrostatic pressure dependence of the Curie 
temperatures. The proposed theory reproduces the experimentally observed linear increase of both 
Curie temperatures with pressure at its low values, as well as the increase of dTci/dp^ at higher 
pressures. The calculated slopes are 3.7 K/kbar at low pressures and 4.3 K/kbar above 12 kbar. 




250 4 ^ 1 ^ > ■ > ■ > 

5 10 15 20 

Figure 7. Hydrostatic pressure dependence of the Curie temperatures of Rochelle salt. Lines: a 
theory; open and closed symbols: experimental points taken from [j^ l and [23 |. 
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To describe the pressure variation of the static permittivity we need to make some changes in 
the theory parameters. The above given value of /ii provides a roughly equal fit to many different 
experimental data for Xii in paraelectric phases at atmospheric pressure. However, the sample to 
sample variation of permittivity at ambient pressure reaches 10% even in good samples [42:j and is 
much larger in samples with defects. This is comparable with the changes in the Curie constants 
produced by hydrostatic pressure [^H^ below 10 kbar. Thus, it seems impossible to try to describe 
these fine high-pressure effects, using for a non-deformed crystal the averaged value of fii given 
in table [T] Instead, we shall determine separate values of ^i, which provide the best fit to the 
experimental data for static permittivity obtained in [2^ and in j4l| at each pressure considered 
therein. Thus, the possible pressure variation of the effective dipole moment ^1 can be inferred. 

The measurements, reported both in [23] and in |4]]|, revealed a decrease of the peak value of 
permittivity at a lower Curie temperature with increasing pressure. This is attributed (2^ to a 
partial clamping of the samples due to the increased viscosity of the pressure-transmitting fluid 
and suppression of the piezoelectric shear strain £4. To calculate the dielectric permittivity of a 
partially clamped crystal we assume that the clamping caused by the viscous fluid under hydrostatic 
pressure is uniform throughout the crystal sample. We assume that, at least, the part of shear strain 
£4 induced by the measuring electric fleld (normally at 1 kHz) is smaller than the one given by 
equation The total strain is then equal to 

e, = e,, + k(^E,-^^,). (23) 



\ ^44 1^*^44 

Here £^4 is a spontaneous part of the strain; is the fleld-induced part of the order parameter. 

The introduced phenomenological coeflicient < fc < 1 describes the extent to which the 
external pressure suppresses the shear strain £4: the cases k — and A; = 1 correspond to the 
totally clamped and free crystals, respectively. Values of k are, naturally, pressure and temperature 
dependent and are different for experimental setups with different pressure-transmitting liquids. 
Thus, the level of clamping at Tci was apparently much higher in the experimental setup of 
(4H than in [2^, possibly because the corresponding temperatures are lower, and the viscosity of 



pressure-transmitting liquid is higher (benzine and silicone oil 4l| vs pentane isopentane mixture 

Substituting equation (|23p into equation (jlOp we obtain the corresponding polarization. Dif- 
ferentiating it with respect to Ei, taking into account equation ([7]) and neglecting variation of 
diagonal strains with the fleld, we get the dielectric susceptibility of a partially uniformly clamped 
crystal 

X^.-X??+ y ^f, , (24) 
2veo (p2 - kAip3 



where 



Xii = Xli + -eUu , /^^ = Ml - 2fc^4d?4 
£0 



From equation in the limiting cases fc = and fc = 1 we obtain susceptibilities of totally 
clamped and free (in the paraelectric phases) crystals. In the ferroelectric phase, a more accurate 
expression for susceptibility should also contain terms like ^eudu, albeit small, produced by 
contributions of diagonal strains. In this subsection these contributions will be neglected, and 
permittivity will be calculated using equation (I24p . 

The temperature curve of dielectric permittivity near the lower Curie point at atmospheric 
pressure presented in (23| appears to be drawn qualitatively. In the fltting procedure we relied on 
the data of Q for these temperatures. Near the upper Curie point, the data of [23| and Q agree 
fairly well. 

Comparison of the calculated temperature dependences of the permittivity with experimental 
data is given in figures [5] and [HI A fairly good description of the experiment in paraelectric phases 
is obtained at a proper choice of k and /ii values. The disagreement at 2.2 and 4.1 kbar in figure [5] 
is due to the mismatch between the calculated and experimental Curie temperatures at these 
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pressures. The disagreement observed in the ferroelectric phase is due to the essential domain 
contributions to the permittivity, not included into the present model. 

To fit the data of the coefficient k is taken 
to be 1 at 297 K (a free crystal). At 255 K we use 
the following values of k: 0.65, 0.4, 0.25, 0.2 at 0, 
0.5, 1.2, 2.0 kbar, respectively. A linear interpola- 
tion between these values at 255 K and 1 at 297 K 
is used. For 3.2 kbar we use fc = at 265 K and 1 
at 297 K, also with a linear interpolation for the 
intermediate temperatures. For the data of [231 we 
use k = 0.9 above 11.7 kbar for the lower Curie 
temperature. 

The pressure and temperature variation of the 
dipole moment /ii is more complicated and contra- 
dictory. First, to fit the experimental data of [4l!] 
for permittivity at ambient pressure we have to as- 
sume that /ii increases with temperature; whereas 
for the data of j23|, Q it should be assumed to 
decrease with increasing temperature, but slower 
than it was given in table [T] 

To fit the permittivity points at high pressures 
we need to assume that the pressure dependence 

of fj,i is opposite to its temperature dependence: if it decreases with increasing temperature, then it 
increases with pressure and vice versa. Thus, for the data of [2^ we use the following dependence 
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Figure 8. The temperature dependence of the 
inverse static permittivity of Rochelle salt 
at different values of hydrostatic pressure ph 
(kbar): 1, □: 0; 2, Q: 0.5; 3, A: 1.2; 4, V: 
2; 5, [>: 3.2. Lines: a theory; symbols: exper- 
imental points taken from [41| . 



Mi=Ai? [l + kr (T-T^,)] (1 + fcpp) 

with = 8.94 • 10^30 (j,^ .^^^ ^ -0.001 K^^. The pressure coefficient fcp was 0.03 kbar^ 
pressures below 5 kbar and above 10 kbar. 
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r(K) 
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Figure 9. (Color online) The temperature dependence of the inverse static permittivity of 
Rochelle salt near the upper (left) and lower (right) Curie points. The values of hydrostatic 
pressure are: 1, *, □: 0; 2, A: 2.2; 3, Q- 3; 4, V : 4.1; 5^ >: 5; 6, o: 11.6; 7, <l: 14; 8, ©: 16.6; 
9, ffl: 20. Symbols are experimental points taken from Q {'k) and [2^ (other symbols). Lines: 
a theory. 



This behavior is quite unusual since the dipole moments are expected to be reduced by hydro- 
static pressure or by decreasing temperature due to the overall reduction of interatomic distances. 
The increase of /xi with hydrostatic pressure can be explained with the hel p of a certain assump- 
tion used in constructing the spatial four-sublattice model of Rochelle salt [4J|. It states that the 
dipole moments in it are actually the 3D vectors, which are not oriented along the a-axis like in 
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the two-sublattice model. Their projections on the b and c axes are different from zero, but com- 
pensated at all temperatures, unless an electric field perpendicular to spontaneous polarization is 
applied. The a-projections are the dipole moments of a two-sublattice Mitsui model, compensated 
in paraelectric phases. It may be assumed that hydrostatic pressure rotates the spatial dipoles 
in such a way that their projections on the a-axis (and /ii) increase. This increase is fast at low 
pressures and slows down above 10 kbar, possibly because the dipoles are already oriented along 
the a-axis. 

To fit the data of [411] we use = 1.0 ■ 10"^" C-m, fcx = 0.003 K-\ and fcp = -0.012 kbar^i. 
This is consistent with the picture when pressure and temperature affect the dipole moment fii in 
the opposite ways. The measurements [41] of humidity (extreme drying and wetting) effect on the 
permittivity of Rochelle salt have shown that in strongly wet samples the permittivity increases 
significantly as compared to the permittivity of normal samples, especially in the upper paraelectric 



phase. It appears that the samples used in the hydrostatic pressure studies 41 ]were moderately wet. 
On the other hand, a decrease of ^ii with increasing temperature (like in [s, 23, 4^]) seems to be an 
intrinsic behavior of normal Rochelle salt samples. To explain the increase of fj,i with temperature 
in wet samples we can assume that due to the excess of water, the crystal conductivity increases, 
thus increasing the dielectric permittivity, especially at high temperatures where the mobility of 
charge carriers is very high. 

4.5. Uniaxial pressures effects 

To get a clearer picture of the uniaxial pressure effect on the Curie temperature and dielectric 
permittivity, it is useful to analyze it within the phenomenological approach. We start with the 
thermodynamic potential 



1 11 



Gi 

V ' 2 ^ -j-'^j 2 

ij—l i—l i—1 



(where sf^ = Mu^Pi are the elastic compliances, and the Mii4 are simply the proportionality 
coefficients between and polarization. Also, Qu are the electrostriction constants; a — aT{T — 
Tq2) or a = aT(?ci ~ '^)' coefScients of the Landau expansion). Hence, one gets the 

following expressions for the temperature and magnitude of the permittivity maxima 

X 2/3 

2 3 I ^ 

= T^l 2 ± — E Q^^P^ + 7 -^^ - E Mu4P^Pi + 514P4 

ax — 4 ut \ ^ — ' 

4=1 y lj=l 

(3 \ 
El-Y Ml^4P^Pi + 914P4 1 • (25) 

In the absence of the fields conjugate to the order parameter {Ei and ^4), the T^^^ correspond 
to the phase transition temperatures. As one can see, in this approximation the uniaxial and 
hydrostatic (ph — Pi = P2 = P3) pressures lead to linear shifts of the Curie temperatures, where 
the permittivity still diverges. In combination with the shear pressure p4 they smear the transitions, 
shift the permittivity maxima, and lower the peaks height. It also follows from equation (j25p that 
if P4 = 



dTck _ Y dTck 
dph (^Pj 



In table |3] we summarize the experimental data on the uniaxial pressure slopes of the transi- 
tion temperatures of Rochelle salt and present our results obtained within the herein developed 
modification of the Mitsui model. The data of [2^ are estimated from the presented therein eii(T) 
curves, each being measured at a single value of pi near the upper Curie temperature. 
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Table 3. Uniaxial pressure derivatives of the transition temperatures of Rochelle salt (in K/kbar). 





(24] 


[25] 


[26] 


this work 


dTci/dpi 


-29 






-26.6 


dTci/dp2 


15 






14.8 


dTci/dp3 


17 






18.0 


dTci/dip2+P3) 






30 ±2 


34.0 


dTc2/dpi 


35 


35.6 




32.9 


dTc2/dp2 


-16 


-18 




-16.5 


dTc2/dp3 


-8 


-6.7 




-8.5 


dTc2/d{p2 + Pa) 






-22 ± 1 


-26.2 



The theory and experiment for the uniaxial pressures agree within 10%, which is close to the 
experimental error [23] . The theoretical dependence of the Curie temperatures on the biaxial pres- 
sure P2+P3 is a little stronger than the experimental one. Also it can be noticed that equation (|26p 
is not fulfilled. 

Independent measurements (25, 26, 41] of dielectric susceptibility of Rochelle salt under different 



uniaxial and biaxial pressures revealed a decrease of the peak values of susceptibility at the Curie 
points as well as smearing out of the peaks. These effects are enhanced with increasing pressures. 

A uniform partial clamping of samples by an apparatus creating the uniaxial pressures can 
explain the observed lowering of the peaks, but not the smearing of the transitions. As an intrinsic 
phenomenon, both these effects can be caused by application of an external field conjugate to the 
order parameter: the electric field Ei directed along the axis of spontaneous polarization or the 
shear stress (74, that is, by the field which induces polarization Pi in the paraelectric phases. In an 
ideal experiment, no uniaxial or biaxial pressure applied along the orthorhombic crystallographic 
axes should act in this way, because piezoelectric coefficients associated with these pressures are 
zeros outside the ferroelectric phase. 

It appears that the observed smearing of the anomalies is an artefact caused by experimental 
errors. We can think of the following factors that in a real experiment can lead to the smearing. 

(i) Stress inhomogeneity. Even a weak inhomogeneity of the applied pressure, partial clamp- 
ing/sample end constraints, surface irregularity result in a non-uniform strain distribution 
over the crystal sample; thus, in different parts of the sample the phase transition is shifted 
to different temperatures. The higher is pressure, the larger is the difference between these 
temperatures, and the more diffuse is the transition, exactly as observed in (4l| . 

(ii) Stray shear stress 174, whose effect would be enhanced by its combination with the uniaxial 
pressures pi , as it follows from equation (j25p . The stress (74 can arise even at a slight misori- 
entation of the samples, as a component of the uniaxial loading intended to create p2 or p^ 
pressures. There can also occur built-in local shear stresses 174 caused by sample defects, e.g. 
dislocations. 

There will be too much uncertainty if we try to take the effect of these factors into account in 
the theory. Thus, we shall not attempt to describe the behavior of the static dielectric susceptibility 
in uniaxially stressed Rochelle salt crystals. 



5. Concluding remarks 

We proposed a generalization of the deformable Mitsui model which along with the piezo- 
electric shear strain £4 takes into account the diagonal strains ei, 62, and £3 as well. In contrast to 
the previous attempt [l^ , in order to incorporate the diagonal strains into this model, the thermal 
expansion is consistently taken into account. 
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In the mean field approximation we find polarization and the strains, as well as thermal, elastic, 
and piezoelectric characteristics related to diagonal strains. For the case of Rochelle salt we suggest 
an elaborated fitting procedure and choose the set of the model parameter values, providing as 
good as possible consistent description of all these characteristics, as well as of external hydrostatic 
and uniaxial pressure effects. 

The expression for susceptibility of a partially clamped crystal is derived in order to describe the 
behavior of the observed susceptibility of Rochelle salt under hydrostatic pressure near the lower 
Curie point. By fitting the theoretical curves to the experimental points, the pressure variation 
of the effective dipole moment fii is estimated. An increase of /ii with hydrostatic pressure at 
low pressures and a decrease with increasing temperature seem to be an intrinsic behavior for 
this crystal. The interaction constants and the asymmetry parameter were found to increase with 
hydrostatic pressure too. This is consistent with the picture, where the dipole moments in Rochelle 
salt are 3D vectors assuming they rotate under pressure in such a way that their projection 
on the a-axis increases. 

Apart from Rochelle salt, the developed modification of the model with appropriate changes 
can be used for consideration of the pressure effects and thermal expansion in other ferroelectric 
crystals described by the Mitsui model. The fitting procedure, however, for each such crystal will 
require extensive experimental data; further measurements will, therefore, be necessary. 

The presented model and the found values of its parameters are a good starting point for 
developing a model description of ferroelectricity in nanosize inclusions of compounds, to which 
the Mitsui model is applicable, grown in a porous matrix [ll, 13, 14 1. 
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Mofle/ib MiL4yT 3 A<3>~OHajibHi/iMi/i p,edpopMau,if\Mv\: o6'eflHaHi/iM oni/ic 
Bn/iMBy aoBHiuiHix tmckIb i Ten/iOBoro posujMpeHHJi b cer^HeTOBiR 
co/ii NaKC4H406 • 4H2O 

A.n. Momfl p.p. J] 

eBi/mbKurfS, I. P. SaneiP 

iHCTMXyT ctiisHKW KOHfleHCOBa HMX CMCTeM HAH yKpaiHI/l, 79011 JlbsiB, By/1. CBEHLliLlbKOrO, 1 

HaL(ioHa/ibHww yHisepcwTeT "/IbBiBCbKa ncnijexHiKa", 79013 JlbBiB, Byn. C. BaHflepw, 12 

SanponoHOBaHO MOfli/i4>iKaL\iKD fle4>opMiBHoV flBoniflrpaxKOBoT MOfle.ni Mii4yV po6iT [Levitskii R.R. et al, Phys. 
Rev. B. 2003, 67, 174112] ra [Levitskii R.R. et al., Condens. IVlatter Phys., 2005, 8, 881], aKa noc/iiflOBHO 
BpaxoBye fliaroHa/ibHi KOMnoHeHTi/i TeH3opa flec)DopMai4ii?i, L140 Bi/iHUKaKiTb nifl pjeto BOBHiiuHix TWCKiB hh BHa- 
cjiiflOK TenjioBoro poBLUnpeHHa. PospaxoBaHO noB'iq3aHi 3 L41/1MH flecfjopiviaLiiaMU renjioBi, n'esoejieKTpi/iHHi Ta 
npy>KHi xapaKTepucTWKW CMcreMH. Bi/iKopi/icTOByKDHi/i sanponoHOBaHy cxeiviy, fl/ia KpucxajiiB ceTHeTOBoT cojii 
3HaMfleHO TaKMM Ha5ip napaiuiexpiB xeopiT, mo 3a6e3neHye saflOBi/ibHe y3rofl>KeHH5i 3 eKcneprnvieHTa/ibHi/iMH 
flaHMMW fljia 3ajie>KHOCTew Tewneparyp Krapi Bifl riflpocxaTi/iHHoro xa oflHOBicHwx TUCKiB, a TaKO>K rewne- 
paxypHMX 3a/ie>KHOCTeM renjioBwx fle4>opMaL4iH, /liHiwHiix Koe4>iL(ieHTiB renjioBoro po3LUMpeHH5i, npy>KHHX 
CTa/iMX cfj i c^, n'e3oe/ieKTpMHHMX KoecjjiLiieHTiB dii i gn (i = 1,2,3). 3ajie>KHOCTi flie/ieKxpwHHoT npoHM- 
KHOcxi Bifl riflpocTaxwHHoro THCKy oni/icaHO 3a flonoMoroio orpHMaHoro b poSori Bnpasy fljia npoHi/iKHOcri 
HacTKOBO saTMCHyxoro Kpwcxajiy. Bi/iJiBJieHO, mo fli/ino/ibHi MOMeHxw xa napawexp acwMexpiT b cerHexoBiCi 
com spocxaKDXb 3 riflpocxaxuHHUM xhckom. 

KyiiOHOBi c/ioea: cefneroBa cinb, Tennoee posini/ipeHHa, riflpocTaTi/iHHi/iiH ti/ick, offHOBicHnCi tuck, MOfle/ib 
M'myi 
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